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Influence of the sign of the refractive index in the reflectivity of 
a metamaterial surface with localized roughness 



o 

(N 
D 

m 



C/3 

o 

• 1— I 

> . 

q 
o 

Oh' 



Vivian Griinhut ^ and Ricardo A. Depine ^ 

Grupo de Electromagnetismo Aplicado, Departamento de Fi'sica, Facultad de Ciencias Exactas y Naturales, Universidad de 
Buenos Aires, Ciudad Universitaria Pabellon I, C1428EHA, Buenos Aires, Argentina 

Received: date / Revised version: January 12, 2013 

Abstract. To study tiie scattering properties of metamaterials, we generalize two scattering methods de- 
veloped for conventional (non-magnetic) isotropic materials to the case of materials with arbitrary values 
(positive or negative) of magnetic permeability and electric permittivity. The generalized methods are used 
to study the changes in the reflectivity of a metamaterial surface with localized roughness when the relative 
refractive index changes sign. Our results show that, unlike the case of a plane surface whose reflectivity 
is unaffected by the change of sign of the relative refractive index, in rough surfaces the change of sign is 
manifested in the reflectivity, even for very low roughness, particularly in observation directions away from 
the specular direction. 

PACS. 42. 25. Fx Diffraction and scattering - 42.68.Mj Scattering, polarization - 42.70.-a Optical materials 
- 78.67.Pt Multilayers, superlattices, photonic structures, metamaterials - 78.68.-l-m Optical properties 
of surfaces 
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1 Introduction 

The area of metamaterials has grown rapidly during the 
first decade of the twenty-first century. While there is still 
no consensus on the definition of the term metamaterial, a 
short and broad definition could be: an artificial environ- 
ment with electromagnetic properties nonexistent or very 
difficult to find in natural materials [T1[5J|3]. Among these 
properties, the one which has perhaps attracted more at- 
tention from the scientific community is the negative re- 
fractive index [3] . In the ideal case of lossless media, neg- 
ative refractive index occurs when there is a range of fre- 
quencies over which both the electric permittivity e and 
the magnetic permeability /i are simultaneously negative 
[S]. For this reason many authors refer to materials with 
negative index as materials with negative constitutive pa- 
rameters [2]. For the actual case of media with losses, the 
condition of negative refractive index is wider (6j and can 
be written as 

efljA^I < 0, (1) 

where en = Re e, fXR = Re fi and Re indicates the real 
part of a complex quantity. It has been recently suggested 
that this condition is only valid for passive media [7]. 

While in a material with positive refractive index the 
vectors electric field E, magnetic field H and direction of 
propagation A; of a plane wave form a right handed set. 
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in a material with negative refractive index these vectors 
form a left handed set. This is why some authors use the 
terms right-handed (RH) and left-handed (LH) to refer 
to materials with positive and negative refractive index 
respectively, even though we think it would be more ap- 
propriate to reserve such designations for chiral media [5] , 
where the internal structure of the medium is associated 
with a direction of rotation. As a result of the positive 
or negative character of the triplet E, H, k, the Foynt- 
ing vector 5' becomes parallel or antiparallel to the wave 
vector k. Therefore, conventional materials are also re- 
ferred to as positive phase velocity (FFV) media whereas 
materials with negative refractive index are referred to as 
negative phase velocity (NFV) media [HIH]- 

The seemingly simple change of sign of the refrac- 
tive index produces dramatic changes in well-known phe- 
nomena as the Doppler effect, the law of refraction or 
Cerenkov radiation [5]. In this paper we are interested 
in investigating the changes in the scattering properties 
of a non-periodic rough surface that separates a conven- 
tional medium from a metamaterial when the refractive 
index of the metamaterial changes sign. Similar problems 
have already been discussed in the following cases: i) lim- 
ited volumes of simple shape, such as cylinders pUlIll) . 
spheres P^fTB] and ii) gratings formed by isotropic [HI 
[T?1ITB1[T7] and uniaxial media [ISiUQj- However, to our 
best knowledge there are no similar studies so far for the 
paradigmatic case of two isotropic half-spaces separated 
by a non-periodic rough surface. Studies of this type may 
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be relevant not only in novel applications where the prop- 
erties of metamaterials play a crucial role, for example in 
the design of invisibility cloaks [^[^ , perfect lenses , 
control of the Casimir force [33] or excitation of surface 
polaritons [^125) . but also in more conventional applica- 
tions, similar to those used for non-magnetic media, such 
as determining the constitutive parameters of a metama- 
terial from experimental reflectance curves as a function 
of the incidence angle [53] . Besides the interest motivated 
by the mentioned applications, investigating the influence 
of roughness on the electromagnetic response of almost 
flat surfaces that differ only in the sign of the refractive 
index on both sides of the surface may also reveal features 
that could be used in nondestructive analysis techniques to 
distinguish the PPV and NPV character of a given rough 
surface. This is so because for perfectly flat surfaces be- 
tween lossless media, reflectance curves as a function of 
incidence angle 6*0 do not distinguish between PPV and 
NPV media with the same absolute value of the relative 
refractive index. This property is part of a more general 
conjugation symmetry [27], valid even for lossy refracting 
media, which ensures that the transformation 

{e^-e\^i^-^i*], (2) 

(where e and ^ now represent the relative parameters and 
the asterisk denotes the complex conjugate) changes the 
phase but not the magnitude of the Fresnel coefficients for 
the amplitudes of reflected and transmitted fields. This 
conjugation symmetry is only valid for non-evanescent in- 
cident waves, i.e, for real angles of incidence and with 
^ 1^0 1 ^ 'k/^. Taking into account that the presence of 
roughness introduces non specular components in the to- 
tal fields generated by an illuminated boundary, the sym- 
metry mentioned before suggests that a far field indicator 
of the PPV or NPV character of a rough surface could 
be found in observation directions away from the specular 
direction. Alternatively, because the presence of rough- 
ness introduces evanescent (non radiative) components in 
the total fields generated by an illuminated boundary, this 
conjugation symmetry, together with the fact that evanes- 
cent waves behave in opposite ways in PPV and in NPV 
media |5,i22!, indicates that the PPV or NPV character 
of a rough surface should also be revealed in near field 
observations. 

In order to explore these issues in an electromagnet- 
ically rigorous way, numerical treatments are inevitable 
and it is convenient and desirable to develop efficient and 
simple methods to evidence the physical mechanisms in- 
volved in the interaction between the incident wave, the 
rough surface and the properties of the refracting mate- 
rial, while preventing the physical mechanisms from being 
masked by the numerical treatments. Since it is not gener- 
ally easy to combine both features in multi-purpose meth- 
ods, such as finite element or finite-difference time-domain 
methods where Maxwell's equations are discretized from 
the beginning, in this work we present two relatively sim- 
ple treatments that meet the above conditions and which 
are based on what is known as the Rayleigh hypothesis 
[2H]. This hypothesis, used by Lord Rayleigh in 1907 to 



solve the dispersion of an acoustic wave by an impenetra- 
ble periodic surface, states that the fields near the cor- 
rugation are composed of waves moving away from the 
surface, an assumption that was objected by Lippmann 
[25] in 1953. Lippmann's objection led to numerous stud- 
ies devoted to establishing the limit of the validity of the 
Rayleigh hypothesis (for an historical overview see for ex- 
ample [30])). The application of this hypothesis in vari- 
ous configurations continues to arouse interest, as shown 
in references [3ni32l[33ll34ll35j . Nowadays, it is recognized 
that in the case of conventional materials the Rayleigh 
hypothesis gives good results for low roughness surfaces. 
The same has been verified in the case of metamateri- 
als with negative refractive index, where the Rayleigh hy- 
pothesis has been successfully employed to study diffrac- 
tion from periodically corrugated surfaces [H1[T51[TB1[T5] . 
In this paper we use the Rayleigh hypothesis to investigate 
the changes in the scattering properties of a non-periodic 
rough surface that separates a conventional medium from 
a metamaterial, when the refractive index of the metama- 
terial changes sign. In Section[2]we provide a brief descrip- 
tion of the boundary value problem for the scattering of a 
plane wave at a rough metamaterial surface, obtaining a 
system of coupled integral equations for the amplitudes of 
the scattered fields on both sides of the surface. Next, we 
decouple the system of integral equations in order to ob- 
tain the generalization for magnetic media of the reduced 
Rayleigh equations (previously obtained by Toigo et al. 
[3 6) for conventional materials) and outline two methods 
of resolution that allow us to calculate the reflected and 
transmitted fields in an independent way: a direct numer- 
ical method, limited only by the validity of the Rayleigh 
hypothesis, and a perturbative method, valid when the 
height of the corrugation is small compared to the wave- 
length of the incident radiation. The perturbative method, 
originally developed by Rice [37 for impenetrable media, 
leads to a relatively simple numerical treatment and is 
very useful to validate the results obtained with the di- 
rect numerical method. Section |3] is devoted to discussing 
the numerical results obtained with both methods for the 
case of deterministic surfaces with a single corrugation, 
postponing the study of statistically characterized rough 
surfaces for future work. We use a time dependence of the 
type e^'"* where uj is the angular frequency, t the time 
and i = \/--T. 



2 Analysis 

2.1 The boundary value problem 

Consider a rough surface represented by the function y = 
g{x) (see Fig.[T]). This surface separates two homogeneous 
and isotropic materials characterized by the constitutive 
parameters (electric permittivity) and /i^ (magnetic per- 
meability), i — 1,2. Medium 1 (y > g{x), medium of inci- 
dence) is a conventional material with positive refractive 
index vi = y/e\Jii, ei > 0, /xi > 0, while medium 2 (y < 
g{x), refractive media) is a metamaterial with frequency- 
dependent constitutive parameters £2 = e2R + 1^21 and 
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Fig. 1. Outhne of the boundary value problem. 



fJ'2 — ^J■2R + *M2/, real parts e2R and i^2R of arbitrary sign 
and positive imaginary parts £2/ > and fi2i > 0. The 
metamaterial, with refractive index 1/2 = ^e2/i2, is either 
an NPV material, if its constitutive parameters satisfy 
eq. ((!]), or a PPV material othervifise. The rough surface 
is illuminated by an electromagnetic, linearly polarized 
plane wave that propagates in the (a;, y) plane (incidence 
plane) and forms an angle Oq, {\9o\ < 7r/2) with the y axis. 
We analyze two independent polarization cases separately: 
the s or TE polarization (electric field in the z direction) 
and the p or TM polarization (magnetic field in the z 
direction). In both cases, the scattered fields (reflected 
and transmitted) conserve the polarization of the incident 
wave. We denote by W(x,y) the z-directed component ei- 
ther of the total electric field (s polarization) or the total 
magnetic field {p polarization). It is known [38 that out- 
side the corrugated region (ming(x) < y < majig{x)), 
^{x, y) can be rigorously represented by superpositions of 
plane waves. If y > u].axg{x), 

-1 /' + 00 

+ ^ / i?(a)e<""+^° '^)da (3) 

represents the incident plane wave (first term, with unit 
amplitude) and the scattered fields in medium I (second 
term, reflected field), while if y < miiig{x), 

-1 /'+00 

•^2(2:, y) = —/ r(a)e<"^-^»'^)da (4) 

represents the scattered fields in medium 2 (transmitted 
field). The quantity ao — kQi^i sin 6*0, ko — w/c, represents 
the X component of the incident wave vector. Note that 
the integrand in represents a plane wave with ampli- 
tude R{a) and wave vector k^^'^\a) = ax + Pa ^y, while 
the integrand in (|4]) represents a plane wave with ampli- 
tude T(a) and wave vector fe''^*''(a) — ax — P^a^y. The 
components along y of the wave vectors fc'^'"^ and fc'^*^ 



are 

|3^J)=|3(^^ia)^{kle,^,,-a'y/\ j ^ 1,2, (5) 

and we define (3^ ^ — /3^^^(ao)- Also, note that the quan- 
tities Pa"^ are real or purely imaginary. In the first case, 
which occurs in the so-called radiative zone |a/A:o| < vi, 
we must require fulfillment of the condition Re ^ > 0, in 
order that the fields in eq. ^ represent propagating plane 
waves that move away from the surface into the half-space 
y > g{x). In the second case, which occurs in the so-called 
non radiative zone |a//co| ^ we must require fulfillment 
of the condition Im ^ > 0, in order that these fields rep- 
resent evanescent waves that attenuate for y — > +00. Sim- 

(2) 

liar considerations are valid for the quantities pa in the 

ideal case of completely transparent transmission media 

(lossless, €2 y /^2 real), although it should be noted that 

in this case the choice of the branches of the square root 

function in ([5]) depends on the PPV or NPV character of 

medium 2. On the other hand, in the real case of lossy 

transmission media (Ime2 > 0, Im/i2 > 0), the quantities 
(2) 

Pa are always complex with a nonzero imaginary part, 

lm/3a > 0, in order that the fields in eq. (HI attenuate for 

(2) 

y — >■ —00. Note that the condition Im/3a ' > automat- 

(2) 

ically sets the sign of Re /3q , independently of the signs 
of e2R and of 112R, i-e., independently of the PPV or NPV 
character of medium 2. 

In order to obtain the unknown amplitudes R{a) and 
T{a) the appropriate boundary conditions at y = g{x) 
must be imposed. These conditions can be written in the 
following way 

>Fi{x,g{x)) ^'F2{x,g{x)), (6) 
j-^^^,ix,gix))^j-^W2ix,gix)), (7) 

where aj — iij for the s mode or aj — ej for the p mode 
and h is the unit vector normal to the surface. 



2.2 Rayleigh-hypothesis 

In order to satisfy the boundary conditions © and ^ 
we use the Rayleigh [28| hypothesis, i.e., we assume that 
the equations ^ and dH), that strictly represent the fields 
outside the corrugated zone imng{x) < ?/ < max(7(x), 
can also be used to represent the fields near the surface. 
Proceeding in this way and projecting the boundary con- 
ditions in the base of Rayleigh functions {e*" ^}a'eRc, we 
obtain a system of two coupled integral equations, whose 
unknowns are the complex amplitudes R{a) and T{a). 
This system can be decoupled through a procedure simi- 
lar to those presented in Refs. [36,39], thus obtaining one 
integral equation for the unknown amplitudes R{a) 

r+00 

-K{ao,a)= K''{a' ,a)R{a')da' (8) 
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and another integral equation for the unknown amplitudes 
T{a) 

- 2/3(1) ^ ^ _ ^ / ^t(^/^ a)Tia')da' (9) 

where 6 () is the Dirac delta distribution. The equations 
(O and (O are Fredholm integral equations of the first 
kind, with kernels 



and 



where 



and 



1 r+°° 

D\u, v\ = — / dxe 



■iux ^ — ivg{x) 



(10) 



(11) 



(12) 



(13) 



is the Fourier transform of e The inhomogeneity in 

eq. is given by K(piQ,a), with 



where 



N 



_ (l-g^)(aoa-/3i^'/3^^')+fcg(ga.^2-.^2) 



delta distributions, we expect the functions R{a) and T{a) 
in the case of slightly corrugated surfaces to be highly 
concentrated around specular observation directions, that 
is a ~ aQ. Neither this feature nor the existence of a 
Dirac delta in eq. © are convenient from a numerical 
point of view, but this can be overcome by introducing 
new fimctions R{a) and T{a) defined as 



R{a) = {a - ao) + R(a) , 

r(a) = r(o),5( 

with and r(o) 



(2) 



T2 qW 



(2) ' 



0(1) , 0(2) ' 



(16) 
(17) 

(18) 
(19) 



the Fresnel coefficients for a perfectly flat surface. The new 
integral equations for the complex amplitudes R(a) and 
f (a) are 



(14) and 



(15) 



K{ao,a) = R^^\a)K''{aQ,a) 

/+00 
K''{a\a)R{a')da', 
-00 



2/3«^=T(«)(a)K*(a,ao) 

p + QO 

+ / K\a,a')f{a')da' . 



(20) 



(21) 



2.3 Direct numerical method 



2.4 Perturbative method 



To solve integral equations as (O numerically, we first 
use a quadrature scheme that allows us to approximate 
the integral as a linear combination of the values R{aj) 
of the unknown function R evaluated at the points of a 
grid {aj}^^i of the independent variable a. The numer- 
ical parameter Na controls the grid density and will be 
determined by convergence criteria. Second, the approxi- 
mated version of the identity ([S]) is evaluated in the dis- 
crete points {aj}fj^i, thus obtaining Na algebraic equa- 
tions whose inversion will allow, in principle, to deter- 
mine the unknowns R{aj), j = 1, . . . , Na- The fact that 
the integration interval of the a variable in is infi- 
nite can be overcome by assuming that |-R(a)| — > when 
|a| — > ±00. In this case, the integral over the infinite in- 
terval can be approximated by an integral over a finite 

interval \a\ < a^^'^, where a^^''^ is another numerical 
parameter to be determined a posteriori through conver- 
gence criteria. A similar treatment for the equation ^ 
allows us to determine T{aj), the values of the unknown 
function T evaluated at the points of the grid {ajjj^i- 

Taking into account that for a flat surface {g{x) = 0) 
the functions R(a) and T(pt) are proportional to Dirac 



When the height of the corrugation is small compared to 
the incident wavelength A, equations (|5]) and ^ can be 
solved by means of a standard [5S1I57] perturbative ap- 
proach. To do so, we introduce the following power se- 
ries expansions for R{a), T(a) and for terms of the form 
^-ivg(x) ^ff\yy2h appear in the kernels pUj) and ([TT|) and in 
the inhomogeneity dT?l) 



R{a) = 



R'^'^\a) 



ri=0 



T(")(a) 



n=0 



-ivg{x) 



E 



(-iwg(a;))" 



The integral can be written as 



(22) 



(23) 



(24) 



(25) 
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where g^'^\u) is the Fourier transform of the function 
[^(a:)]" and the index n in the series ([22]), and ([25]) 
indicates the perturbative order. When these expansions 
are introduced in the integral equations (O and (O the fol- 
lowing iterative schemes for the coefficients i?'^"-'(a) and 
T(")(a), n > 1, are obtained 



(26) 



a 



En 



^(J)(Q,_Q,')T("-J)(a') 



where R'^'^\q.) and r*^°^(a) coincide with the Fresnel coef- 
ficients given by the equations (fTSi) and (|T9| . 



3 Numerical treatment 

The methods presented in the previous section have been 
implemented numerically for surfaces with a finite number 
of protuberances limited to region —a/2 < a; < a/2. To 
illustrate the changes produced in the reflectivity of a non 
flat surface when the sign of the relative index of refraction 
is changed, the incidence medium is vacuum (ei = 1, /ii = 
1) and the transmission medium is a lossy metamaterial in 
all examples presented here. The values of the constitutive 
parameters of the metamaterial are £2 = 5-1- 0.01 i, ^2 — 
1-1-0.01 i (a PPV medium with refractive index V2 ~ 2.23-1- 
0.01 i) or es = -5 -h 0.01 /i2 = -1 + O.OH (an NPV 
medium with refractive index V2 ~ —2.23 -I- 0.01 1). Bear 
in mind that one set of constitutive parameters is obtained 
from the other set through the conjugation transformation 
expressed in eq. ^ and therefore both sets give the same 
reflectivity for a perfectly flat surface. 

To control the quality of the calculations, we have 
checked the convergence of the results for different val- 
ues of the numerical parameters and the agreement be- 
tween the direct and the perturbative methods. Besides, 
as energy is conserved in the scattering process, we have 
checked the fulfillment of the power conservation criterion. 
Taking into account that we are considering lossy media, 
it is convenient to write this criterion in the following form 
001 

+ = 1 , (28) 



where 



Pr = 



Rc 
2^ 



+00 ^(1) 



(lyl^(a)l da, 



(29) 



represents the fraction of the incident power which is scat- 
tered (reflected) into the incident medium, and 



Pn 



CTl Rc f 
-2 J- 



+ 00 



I dy dx 3 



V=a(x) 



dx , 



(30) 



represents the fraction of the incident power which is ab- 
sorbed by the medium below the surface. 

To evaluate the reliability of the application of Ray- 
leigh methods to non periodic rough surfaces with nega- 
tive refractive index, we consider in this section the case 
of a rectangular protuberance (width a and height /i), il- 
luminated at normal incidence (6*0 = 0°). In this case, 
g{x) = hmcix/a), where rec(u) is the rectangular func- 
tion centered at the origin with unit width and height. 
In Figures [5] (for s polarization) and [3] (for p polariza- 
tion) we compare the curves |i?(Q;)p vs a/feg obtained with 
the direct method (continuous curve) and the first-order 
perturbation method (circles) for the case h/a = 0.005, 
A/a = 0.5 and for PPV and NPV media. We observe 
that both methods give an excellent agreement and that 
in all the cases they predict the presence of a principal 
maximum centered at the value of the spectral variable 
a = QfQ, that corresponds to the specular reflection di- 
rection. This maximum, of width 2A/a, is surrounded by 
secondary maxima of width A/a, in total coincidence with 
the results predicted by the scalar theory of diffraction for 
a slit 2A wide. The power conservation criterion [28] is sat- 
isfied with an absolute error less than 0.001. These results 
show that when the height of the protuberance is small, 
the perturbative method converges rapidly and coincides 
with the direct numerical method, as it has been already 
observed for periodical corrugations in NPV media |14) . 
We have verified that the same occurs for other values of 
geometric and incidence parameters and for other protu- 
berance shapes. 

In Figures [?] (s polarization) and [5] {p polarization) 
we repeat the comparison between the direct (continu- 
ous curve) and the perturbative (circles) methods, for 
the same situations considered in Figures [2] and [3] ex- 
cept that now the height of the protuberance is 10 times 
higher {h/a = 0.05). It is interesting to observe that in 
this case the perturbative method converges better in the 
PPV case, where the coincidence with the direct method 
is obtained in both polarizations in the 8*'' perturbative 
order, while in the NPV case the coincidence is obtained in 
the 20*'* perturbative order. For this protuberance height 
the power conservation criterion[5S]is satisfied with an ab- 
solute error less than 0.07 (PPV) or less than 0.13 (NPV). 
We have found that for greater heights and for NPV media 
the perturbative method can fail to converge, while it still 
converges when the medium is PPV. For example, for a si- 
nusoidal protuberance of the form g{x) — ^[1 -|-cos(^a;)] 
rec(a;/a), h/a = 0.1, A/a = 0.5, 6*0 = 0° and PPV media, 
the perturbative results converge to those obtained with 
the direct numerical method in the 10*'' order, although 
convergence is not obtained when the medium is NPV. 
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Fig. 2. (Color onhne) Curves of |_R(a)p vs a/ko obtained 
with the direct method (continuous curve) and the first order 
perturbative method (circles) for a rectangular protuberance 
of width a and height h/a — 0.005 illuminated at normal inci- 
dence {6o — 0°) in s polarization. The wavelength is A/a = 0.5. 
The media are NPV (above) or PPV (below) with the same 
module of the refraction index. 

4 Results 

Having checked the reliability of the results obtained with 
the Rayleigh methods, we now turn our attention to the 
reflectivity of metamaterial protuberances whose indices 
of refraction have opposite signs. The normalized angular 
distribution of power scattered into medium 1 is given by 
the integrand in eq. ([29| 

^ = -^Re^|i?(a)p, (31) 
da zna p^'-' 

an expression which shows that only the values of R{a) 
in the radiative zone, |a/fco| < I'l, contribute to the scat- 
tered power. In this spectral zone the integrand in eq. ([3]) 



0,0100-, 




Fig. 3. (Color online) Curves of |_R(a)p vs a/ko obtained 
with the direct method (continuous curve) and the first order 
perturbative method (circles) for a rectangular protuberance 
of width a and height h/a = 0.005 illuminated at normal inci- 
dence {00 = 0°) in p polarization. The wavelength is X/a — 0.5. 
The media are NPV (above) or PPV (below) with the same 
module of the refraction index. 

represents plane waves propagating away from the sur- 
face along a direction that forms a scattering angle 9 si, 
{\Osi \ < 7r/2), with the +y axis. 

In Figure [H] we compare the angular distributions of 
reflected power that correspond to a single sinusoidal pro- 
tuberance illuminated under normal incidence with X/a = 
0.25 and h/a — 0.0025. Wc have verified that the power 
conservation criterion is satisfied to an error less than 
0.2%. As in the examples considered in Figures [5] and [31 
close to the specular direction the four curves (s PPV, 
s NPV, p PPV, p NPV) show essentially the same re- 
sponse and so we have chosen an scale that amplifies the 
differences between the curves in observation directions 
^si 7^ ^0- The first perturbative order exhibits analytically 
the coincidence of the responses in the specular direction. 
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Fig. 4. (Color online) Curves of |-R(q:)P vs a/ko obtained 
with the direct method (continuous curve) and the first order 
perturbative method (circles) for a rectangular protuberance of 
width a and height h/a — 0.05 illuminated at normal incidence 
{9o = 0°) in s polarization. The wavelength is A/a = 0.5. The 
media are NPV (above) or PPV (below) with the same module 
of the refraction index. 
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Fig. 5. (Color online) Curves of |7?(q)P vs a/ko obtained 
width the direct method (continuous curve) and the pertur- 
bative method's first order (circles) for one rectangular protu- 
berance of wide a and height h/a = 0.05 illuminated at nor- 
mal incidence {9o = 0°) in p polarization. The wavelength is 
A/a = 0.5. The media considered are NPV (above) or PPV 
(below) with the same module of the refraction index. 



as can be seen from the expression (j26p . From a physical 
point of view, such coincidence can be understood taking 
into account the low height considered in this example 
and the fact that for a flat surface: i) the s and p po- 
larizations are indistinguishable in normal incidence and 
ii) the exchange between PPV and NPV refractive media 
changes the phase but not the absolute value of the Fres- 
nel coefficients for the amplitudes of the reflected flelds. 
However, it should be noted that despite the low height 
considered, for observation directions away from 9q, 
the angular distribution of the reflected power correspond- 
ing to each incident polarization is sensitive to the sign of 
the refractive index of the metamaterial. Similar features 
are observed for oblique incidences, as shown in Figures [7] 



and [51 obtained for the same parameters used in Figure [SI 
except that now 9q = 20°. 

In Figures [S] (for normal incidence) and [TUj (for Oq = 
20°) we show the results corresponding to a single sinu- 
soidal protuberance with a height value ten times greater 
than the value used in Figures [ni[7| and [5] We observe that 
when the protuberance height increases, the power scat- 
tered near the specular direction {6si ~ 9o) becomes more 
sensitive to the PPV-NPV interchange. 

The results presented in Figures [51 [3 [HI [D and [TUj con- 
flrm the conjecture suggested by the conjugation symme- 
try, that is, the optical response along observation direc- 
tions away from the specular direction could be used as a 
far-fleld indicator of the PPV/NPV character of shallowly 
corrugated metamaterial surfaces. Note that this indicator 
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Fig. 6. (Color onhne) Angular distribution of power scattered 
into medium 1 for a sinusoidal protuberance. The incidence 
parameters correspond to those of Figures [5] and [3] and the 
geometric parameters are A/a = 0.25, h/a — 0.0025. 
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Fig. 8. (Color online) Same as in FigureO but with a difi^erent 
scale to evidence the differences between curves in observation 
directions far away from the specular direction. 



-30 30 




Fig. 7. (Color online) Angular distribution of power scattered 
into medium 1 for a sinusoidal protuberance. The geometric 
and incident parameters correspond to those of Figure [SI ex- 
cept that now Oq — 20°. 



involves only the radiative range of the quantity R{a), i.e., 
the range in which R(a) represents the amplitude of prop- 
agating plane waves. However, in applications where near 
fields are involved, the non-radiative range of the quantity 
R{a) can also play an important role as a PPV/NPV in- 
dicator. This is due to the fact that in the non-radiative 
range the quantity R(a) represents evanescent waves that 
only affect the value of the fields near the surface and it is 
well known that the behavior of evanescent waves changes 
dramatically depending on the sign of the refractive index 
of a medium. To explore this possibility, in Figure [TT] we 
compare the curves of |i?(Q;)p vs a/fco in the non-radiative 



Fig. 9. (Color online) Same as in Figure [51 except that now 
h/a = 0.025. 



zone for PPV and NPV corrugated metamaterials with 
the same geometric and incident parameters considered 
in Figures O It can be seen that, although the surface is 
optically almost flat {h/X = 0.1), the amplitudes of the 
evanescent fields generated during the scattering process 
are sensitive to the change of sign of the refractive index. 



5 Conclusion 

We have extended two electromagnetic scattering formalisms, 
originally developed for conventional (nonmagnetic) ma- 
terials, to the case of corrugated materials with arbitrary 
(positive or negative) constitutive parameters. We are plan- 



Griinhut and Depine: Influence of the sign of the refractive index in the reflectivity of a metamaterial with roughness 9 



1,2 

dP'"/da 



0,0 



p NPV 




• p PPV 




s NPV 




s PPV 






\k 
















h 'A 















-80 -60 -40 -20 20 40 60 80 



0si C^a) 



Fig. 10. (Color online) Same as in Figure (5] except that now 
6*0 = 20°. 



terials and have centered our attention on the scattered 
fields reflected into the medium of incidence. In a next 
step, we are planning to consider the ideal case of loss- 
less metamaterials in order to study the behavior of the 
scattered fields transmitted into the medium below the 
surface. Starting from a perfectly flat surface whose reflec- 
tivity is unaffected by the PPV/NPV transformation, we 
found that the angular distribution of scattered power cor- 
responding to the same surface with a shallow added cor- 
rugation is more sensitive to the PPV/NPV interchange 
when the observation is not close to the specular direction 
and that the sensitivity increases when the height of the 
corrugations increases. We have also shown that the near 
fleld can be highly affected by the PPV/NPV interchange, 
even for surfaces with very low height protuberances. 

We acknowledge financial support from Consejo Na- 
tional de Investigaciones Cientfficas y Tecnicas (CON- 
ICET), Agenda Nacional de Promocion Cientifica y Tec- 
nologica (BID 1728/OC-AR PICT-11-1785) and Univer- 
sidad de Buenos Aires (UBA). 




a/kg 

Fig. 11. (Color online) Curves of |_R(q)P vs a/ko for a single 
sinusoidal protuberance. The geometric and incident parame- 
ters correspond to those of Figures [9l 



ning to use this extension to study the electromagnetic 
scattering from rough surfaces in applications motivated 
by the recent emergence of metamaterials with negative 
refractive index. We have presented examples showing that, 
even though they require different numerical treatments, 
both formalisms give coincident results for shallow corru- 
gations (the known range of validity of Rayleigh methods) 
and that such results are in agreement with the predictions 
of physical optics. We have used these formalisms to illus- 
trate the changes produced in the scattering properties of 
a surface with a protuberance when only the sign of the 
refractive index of the metamaterial below the surface is 
changed. To be realistic, we have considered lossy metama- 
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